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We consider the use of quantum noise to characterize many-body states of spin systems reahzed 
with ultracold atomic systems. These systems offer a wealth of experimental techniques for realizing 
strongly interacting many-body states in a regime with a large but not macroscropic number of 
atoms where fluctuations of an observable such as the magnetization are discernable compared to 
the mean value. The full distribution function is experimentally relevant and encodes high order 
correlation functions that may distinguish various many-body states. We apply quantum noise 
^ . analysis to the Ising model in a transverse field and find a distinctive even versus odd splitting in 

the distribution function for the transverse magnetization that distinguishes between the ordered, 
critical, and disordered phases. We also discuss experimental issues relevant for applying quantum 
fSJ . noise analysis for general spin systems and the specific results obtained for the Ising model. 

I. INTRODUCTION 
5-H ■ A. Strongly correlated systems of cold atoms 
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A natural regime for ultracold akali gases is when interactions are weak. For typical experimental densities of 
n K, 10^^ atoms per cubic centimeter and scattering lengths of a few nanometers, interparticle distances are two 
orders of magnitude larger than the scattering length. This regime was the focus of most experimental work done 
in the first decade since the observation of Bose-Einstein condensation Q, Q, |^ . Excitement in this field was fueled 
(~| , by the possibility of twcmngQedankenexperiments into reality. This includes interference of independent condensates 
O ' 01 , creation of atom lasers 0, , and condensation of multicomponent systems |E S S IS Till IT^ . The Gross- 
O , Pitaevski equation captures all essential features of systems in this regime 0, 0| . This is a mean- field approximation 
that only requires the knowledge of a single wavefunction describing a macroscopically occupied state describing the 
condensate. Conceptually, the many-body aspect is important only in turning quantum coherence into a macroscopic 
^ phenomenon. However, solutions to the non- linear equations may be non-trivial with the vortex lattice [isl ll^ 
OO being a prime example. 

\l I The current focus for the field of ultracold atoms is shifting towards creation and study of systems with strong 
interactions and correlations. The driving force of this evolution is remarkable experimental progress which provide 
powerful tools for realizing strongly interacting systems of cold atoms. Magnetically tuned Feshbach resonances (see 
|l8l ITgl l and references therein) can be used to make scattering lengths comparable to interparticle distances. This 
' allowed creation of molecules and observation of pairing in fermionic systems. Confining particles in optical lattices is 
another powerful tool for amplifying the role of interactions. Atoms in a periodic potential created by a standing wave 
Ci I laser beam can exhibit both enhanced interactions and suppressed kinetic energy. The relative strength of the two 
energies can be tuned by changing the intensity of the laser beam. This technique was used to observe the superfluid 
to Mott transition in which the behavior of atoms changed from wavelike in the superfluid phase to particle like in 
the Mott phase [23|. Reduced dimensionality is another useful tool for making strongly correlated systems. The 
strength of an effective one-dimensional interaction can be changed by tuning either the particle density, transverse 
^ confinement, or scattering length. This approach was used in several experiments that explored the crossover from 
weakly interacting quasicondensate to the Tonks gas of hard core bosons |^ . There are many other promising 
techinques for obtaining strongly interacting systems of quantum degenerate particles including rotating condensates 
|23l | , ultracold polar molecules (see and references therein) , and Rydberg atoms [2^ . 
^ ■ One of the primary motivations for studying strongly correlated systems of cold atoms is the expectation that they 
will realize important and not yet fully understood quantum many body states. The target list of modern experiments 
includes the Fulde-Ferrel-Larkin-Ovchinnikov phase of unbalanced fermion mixtures with attraction |2^|2^ . supersolid 
phase US 113,1^, quantum magnetic phases of atoms in optical lattices (ordered or spin liquids) 32, 33, 34,'35j, 
fermions with d-wave pairing in lattice systems with repulsive interactions |3a | , and quantum Hall phases .37,. 38, 39|| . 
An important open question is detection and characterization of these states. 

Analysis of correlation functions is a common way of characterizing many-body states in condensed matter physics. 
Neutron and light scattering experiments can be used to measure spin and charge response functions. Spin and charge 
density wave orders give rise to additional 5-function peaks in the static correlation functions, hence they immediately 
show up as new Bragg peaks in elastic neutron scattering experiments. Conductivity measurements at finite frequency 
(microwave, infrared, optical, etc.) provide access to current-current correlation functions (4fl||. They are a useful tool 
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for detecting states with a gap in the quasiparticle spectrum, such as superconducting and various density wave 
states. Photoemission experiments measure electron spectral functions and can be used to study excitation spectrum 
or observe electron fractionalization in one dimensional systems 0, • So what we want for strongly correlated 

systems of cold atoms is to have a tool box for measuring correlation functions. 

Certain experimental techniques used in cold atoms provide natural analogues of condensed matter experiments. 
For example, Bragg spectroscopy is similar to light and neutron scattering in solid state physics and can be used to 
analyze the dynamic response functions 'S]. RF spectroscopy in cold atom systems is analogous to finite frequency 
conductivity measurements ,4^ ,46, 47, . There are also several proposals for designing experiments that mimic 
solid state methods |4^ . Additional techniques such as time of flight and interference experiments 0, 13 or the 
possibility of single atom detection |^ are also available which are not feasible with traditional condensed matter 
systems. The richness of physical phenomena that we expect from strongly correlated systems of cold atoms provide 
strong motivation to think about new methods and approaches that would provide "smoking gun" signatures of 
various many-body phases. 



B. From noise analysis to correlation functions in interacting cold atoms systems 

The usual philosophy of performing measurements in physics is to minimize the noise by averaging over many 
experiments. However, there are notable exceptions. In mesoscopic solid state physics and quantum optics, one is 
often interested in the noise rather than the average signal. For example, shot noise was used to measure fractional 
charge and statistics of excitations in fractional quantum Hall systems |53l l54| . Quantum noise in photon counting 
experiments was used to identify non-classical states of light |55j| . 

The philosophy of noise analysis was implemented in several recent types of experiments with cold atoms. These 
experiments analyzed quantum noise in time of flight images of atoms released from an optical lattice [53.15^ Is^lssHsoj l 
(see Ref. pO| for theoretical analysis). Free expansion during the time of flight maps momentum distribution of atoms 
inside the interacting system into the density distribution after the expansion. Correlation functions {p{'r)p[r')) of 
the images in real space contains information about the {n^nki) correlation function of the interacting system in 
momentum space. Note that jik — c^Cfe is the occupation number of momentum states and not the density operator 
at wavevector k. Hence (rifc^fe') represents a non-local correlation function for the original interacting system. 

What can one learn from {nkUk') correlation functions? The simple answer is the periodiciy of the original system. 
Quasimomentum (i.e. lattice momentum) is a good quantum number in a periodic system. Quasimomentum differs 
from the physical momentum by Bragg reflections by reciprocal lattice vectors. This leads bunching for bosons and 
peaks in (nferifc/) when the two momenta differ by a reciprocal lattice vector k — k' = G. Analogously, fermions 
exhibit antibunching. Such experiments have been done for spinless bosons [s^ Is^ Is^ as well as fermions [s^ 
and provided one of the most elegant realizations of Hanbury-Brown-Twiss experiments with neutral atoms. When 
such experiments are performed with mixtures, we can expect even more exciting results. They should provide an 
ideal probe for antiferromagnetism or static density wave, since both phenomena increase the size of the unit cell 
and give rise to additional reciprocal lattice vectors It is worth pointing out that fluctuating orders, like in 

Luttinger liquids, can also be observed using this method |6lj| . Another beautiful application of the noise analysis in 
time of flight experiments was the measurement of fermion pairing on the BEC side of the Feshbach resonance |58j | . 
These experiments demonstrated correlations in the occupation numbers for fermions with momenta k and — fc, which 
demonstrated the existence of a Cooper pair condensate at zero momentum. 

A different method for analyzing correlation functions in interacting many-body systems is based on interference 
experiments between two independent condensates. It may not be immediately obvious why this measurement cor- 
responds to the noise analysis. In every shot we should find an intereference pattern. However the position of 
interference fringes is unpredictable from shot to shot. If we average over many experiments we will find that the 
interference pattern washes out completely. We can put this in a more formal mathematical language: Let A be the 
quantum mechanical operator that describes the amplitude of interference fringes [s^, we have {A) = but finite 
What is interesting for our purposes is the value of itself, or the related quantity of the fringe contrast. 

When individual condensates are not fluctuating, we expect fringes with perfect contrast. When there are strong 
thermal or quantum fluctuations we expect suppression of the fringe amplitude. Theoretical analysis performed in 
Ref. (621] showed that the amplitude of interference fringes can be related to the instantaneous two point correlation 
function: = fl G^{r), where G{r) = (c(r)c^(O)) and Q. is the imaging area. By analyzing scaling of with 

the size of the images region, Hadzibabic et al. recently demonstrated the Berezinskii-Kosterlitz-Thouless transition 
in two dimensional condensates 63] . 

The common philosophy of HBT experiments with atoms [s^, EE 113, IH, [s^ [63, 0, 112, 113] and analysis of the 
average contrast in interference of independent condensates (HC) experiments is the focus on second order coherence. 
This means that in the analysis of absorption images, we are not looking for peaks in {p{r)) but in {p(r)p(r')) |86j ]. An 
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interesting question to ask is whether one can use quantum noise measurements to analyze higher order coherence. In 
the case of IIC experiments this problem was considered in Refs. [6^l6^ . The basic idea is that we do not only measure 
the average value of the amplitude of interference fringes but also their fluctuations. And it is straightforward to show 
that that higher moments of \A\'^ correspond to higher order correlation functions. In the case of one dimensional 
Bose systems with interactions it is possible to find the explicit form of the distribution function of \A\'^. For weakly 
interacting bosons phase fluctuations within condensates are small and one finds the Gumbel distribution of |v4p with 
a small width. When interactions between atoms become so strong that the system approaches the limit of hard core 
bosons, the distribution function becomes Poissonian. 

C. Quantum magnetism with cold atoms in optical lattices 

The main subject of this paper is to extend the approach of quantum noise analysis to the study of higher order 
correlations in quantum magnetic systems realized using cold atoms. Several approaches have been discussed recently 
for creating magnetic systems with cold atoms (see Ref. 35] for a review). Our proposal is sufficiently general so we 
expect that it can be applicable to all of them. The main idea can be understood in the following example. Consider 
a measurement of the a component of spin magnetization in a one dimensional spin systems 

M„=^^„(i). (1.1) 

i 

The average magnetization can be obtained by averaging over many experiments 

{M^)=J2{Sc.{^)) (1.2) 

i 

but individual measurements will exhibit quantum noise. Each measurement will give a value of Ala that is different 
from the average (Mq,). It is easy to see that the strength of quantum fluctuations is controlled by correlation functions 
of the spin operators. For example 

(Ml) - J2 iSc.in)Sa{^2)) (1.3) 

ii ,i-2 

So the width of the distribution is determined by two point correlation functions. By the same argument higher 
moments provide a measure of higher order correlation functions 

(^4">= E {Sa{ll)...Sa{ln)) (1.4) 

The knowledge of high moments is contained in the distribution function of Mq,. So if we know the distribution 
function of P{Ma) we obtain information about correlation functions of arbitrarily high order. In real experiments, 
of course, technical noise and a finite number of experiments will limit the number of moments that can be reliably 
extracted. 

Generalization of this idea to other measurements is straightforward: one should measures not just the average 
value of some operator but its entire distribution function. The distribution function contains information about 
high moments which can be related to high order correlation functions. This provides valuable characterization of a 
many-body state. It is important that measurements are done on a mesoscopic system, hence fluctuations will not be 
neglibly small compared to the mean value. For macroscopic solid state samples this is not easy to achieve, but for 
cold atoms systems with tens or hundreds of thousands of atoms this is an experimentally relevant regime. Another 
advantage of the cold atoms systems is their nearly complete isolation from the environment. We point out that Ref. 
|65| considered a related idea of measuring the number of particles in a mesoscopic part of a fermionic system with 
pairing. They showed that fluctuation in the total number of particles will exhibit interesting evolution as one tunes 
the system from the BCS to the BEC regime. 

In the main part of this paper contained in Sec. ^we consider a specific example of the magnetization distribution 
for the Ising model in a transverse field. We are not the first ones to address this question. Ref. ^fl] study the 
distribution functions for both the longitudinal and transverse magnetizations of this model with analytical and 
numerical techniques. Wc focus primarily on the transverse magnetization but provide a more detailed analysis of its 
global structure. In particular, we study how higher order moments encode the different correlations in the ordered, 
critical, and disordered phases and give a simple physical picture in terms of confinement versus proliferation of 
domain walls. Sec. IIIII discusses experimental issues involved in applying quantum noise analysis to ultracold atoms 
including those specific to the results we obtained for the transverse Ising model. We summarize our discussion in the 
conclusion of Sec. IIVI The appendices contain mathematical details of the discussion in Sec. 
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II. MAGNETIZATION DISTRIBUTION IN THE TRANSVERSE ISING MODEL 



The one-dimensional Ising model in a transverse field provides a textbook example of a symmetry-breaking quantum 
phase transition |67j . It is described by the Hamiltonian 

N 

H = -JY, [2S^{i}SA^ + 1) + gS^ii)] (2.1) 

where Sa{i) are spin-1/2 operators at site i. Here g is a dimensionless coupling which sets the strength of the transverse 
field and J sets the overall energy scale. We take periodic boundary conditions and N is the number of sites. The 
competition between the transverse field and the Ising term of the Hamiltonian drives a quantum phase transition 
at = 1 from an ordered phase g < Qc with spontaneous magnetization along the x direction to a disordered phase 
g > gc with zero magnetization along x. We are interested in the distribution function for the magnetization Ma 
given by 

P(Af„) - ((5(M„)) (2.2) 

where 5{x) is the delta function and brackets denote ground state expectation values. 

The global structure of P{Mx) describing the order parameter follows from low-order moments of the longitudi- 
nal magnetization Mx and symmetry breaking arguments. Away from the critical point g ^ gc, exact results on 
{Sx{i)Sx{j)) '6^, 'go'I show the spin correlation length and the second moment (AfJ) are finite. Physically, this means 
spins fluctuate independently on length scales larger than the correlation length. Thus, for a sum of a large number 
of spins, the corresponding distribution should be approximately gaussian. In the broken symmetry phase g < gc, the 
distribution should have two peaks describing the formation of spontaneous magnetization for Mx- On the other hand, 
the unbroken symmetry phase g > gc should only have one peak centered at zero. At the critical point g — gc, the 
power-law decay of {Sx{i)Sx{j)) ~ I* ^ gives a divergent contribution to (M^) and a non-trivial distribution 

function. 

More analysis is needed for the transverse magnetization where symmetry breaking arguments do not apply. 
Exact results 68, 70] for {Sz{i)Sz{j)) show the second moment (Af^) is finite even at the critical point. Analogously, 
P{Mz) is also gaussian. However, additional aspects of the global structure such as the number of peaks and other 
features are encoded in higher order moments. 

Instead of analyzing individual moments, we study the entire distribution P{Mz) directly. We find its global 
structure refiects the amount of disorder in the system. This is in contrast to P{Mx) which reflects the amount of 
order. We show P{Mz) is asymptotically gaussian for a large number of spins with a characteristic splitting between 
even and odd magnetization that differs between the phases. A duality mapping shows this behavior has a physical 
interpretation in terms of proliferation, criticality, and confinement of dual domain walls when g < gc, g = gc, and 
g > gc, respectively. The domain wall physics is particularly clear in the small and large g regimes where P{Mz) 
describes the qualitatively different counting statistcs of pairs and single domain walls, respectively. In section Hi Al 
we calculate P{Mz) via fermionization. We obtain and discuss its global structure in section ITlBl focusing on the 
even versus odd splitting and its interpretation using a duality mapping. Appendix ^ discusses the use of Toeplitz 
determinants employed in the calcuation. In addition, appendix ^ derives the small and large g behavior in more 
detail. Finally, Appendix El describes the saddle-point analysis used to obtain the asymptotic gaussian behavior. 

A. P{Mz) Distribution Function 

Jordan- Wigner fermionization (see Ref. [tJ) maps spins in one dimension onto spinless fermions 

= X{[~^.{3)]^-{^ (2.3) 

j<i 

c\ = n[-'^^(*)]'^+(*) (2-4) 



where Uaii) = 2Sa{i) are the Pauli matrices and a± = ax i i<Jy are the corresponding raising and lowering operators. 
The spin Hamiltonian of Eq. l2.1l maDS to free fermions 

-ff = [(g + cosfc)(4cfc -f c^fcC^fe) - isinfc(c_fcCfe - c\c^_k) (2.5) 

fc>o 
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where Ck are in momentum space. 

Instead of studying directly, the analysis is simpler for the quantity 



" / 1 \ " 



which is simply related to both and the fermion number. Recall N is the number of spins in the entire chain 
and compare this to n which is the number of spins in a block embedded in this system. We will take iV — > c» first 
to describe the thermodynamic limit and then n — > c». Physically, we are studying the local magnetization Mz of a 
large block of spins embedded in a much larger system. This has important consequences in the interpretation of the 
results and experimental detection. 

We will calculate the generating function for p„ defined as 

Xni\) = {e^''-) (2.7) 

where the brackets denote ground state expectation values in the thermodynamic limit N ^ oo. The corresponding 
distribution function is given by the fourier transform of the generating function 

fix 

Pn(m) = / —e-'^"^Xn{X) (2.8) 



.. 27r 

which is related to P{Mz) by the following 

P{Mz)^Pn{nl2-Mz). (2.9) 
We first define the following linear combinations of fermionic operators 

A, = c, + c\ (2.10) 

B,{X) = c, + e'^c\ (2.11) 

and obtain for the generating function 

X„(A) = (AiBi(A) . . . A„B„(A)) (2.12) 

where we have used the operator identity e*^(i-c]ci) _ AiBi{X). Using Wick's theorem, the expectation value in Eq. 
12. 121 can be reduced to products of the following two-point functions 

{A:A,) ^ r ^e-^^-(^-^) (2.13) 

^^-^k(^-J)^-2^X ^ g.^(^) ^j^^^g^)] (2.14) 

(ABjiX)) = r ge-*'=(^-^)i [(1 + e'') + (1 - e'^)e'^^] (2.15) 



where the Bogoliubov angle 



+ e*'' 

^ (2.16) 



comes from diagonalizing the Hamiltonian in Eq. 12.51 using a Bogoliubov transformation. 

Notice (AiAj) = Sij. However, the contractions of this form entering Xn(A) vanish because they always occur with 
i 7^ j. Moreover, any full contraction containing {B.iBj) also vanishes since it must contain some (AiAj) contraction. 
This implies the only full contractions that contribute are those that only pair Ai with Bj. Taking into account the 
fermionic signs, these full contractions organize into a determinant of the form 

X„(A) = .det [/.-,(e^^)] (2.17) 

t,j — l...n 



where fn{z) are the Fourier coefficients of the function 



(l + z) + (l-z) 



V(3 + a:)(g + x-i) 



(2.18) 
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Notice that the coefRcients are given by /n(e*'*') = {AiBj{\)) and the function f{z,x) is the integrand in Eq. 12.151 
with z — e^^ and x = e^^ . 

The above resuh gives an exphcit expression for the generating function. Here we can make some general comments 
on the structure of Xn(A). By diagonahzing the matrix inside the determinant, we see that the generating function is 
of the form 

n 

Xn{X)=X[^[[l + e^'') + [l^e^'')^^^]■ (2.19) 
1=1 



My = f ^e-'k{^-3)e^o,^ (2.20) 



Here fii are the eigenvalues of the n x n matrix Mij formed from the Fourier coefficients of the Bogoliubov angle 

r 

„ 27r 

This can be rewritten in the more suggestive form 

n 

Xn(A) - llxB{X,P^) (2.21) 

i=l 

XBiX,p) = (l-p)+pe*^ (2.22) 

Pr = ^(l-M^) (2.23) 

as a product of generating functions xb{X,p)- We recognize these generating functions as that of a Bernoulli random 
variable with probability p of obtaining one and 1 — p of obtaining zero. This implies pn is the sum of independent 
Bernoulli variables. This is not surprising since p„ is directly related to the fermion number and each fermionic level 
can be either occupied or unoccupied. 

However, even though the matrix is real, the eigenvalues /i^ are only required to occur in complex conjugate 
pairs. In fact, we find the fii are generically complex and satisfy < 1. This means that XB{X,pi) cannot be 
interpreted as a classical Bernoulli distribution. However, a pair of conjugate eigenvalues gives 

XB{X,Pr)XB{X,P:) = i [|1 - + 2(1 - \M')e'^ + |1 + (2-24) 

allowing an interpretation as a classical distribution with the coefficient of e™"^ giving the probability of obtaining to. 
Notice that the interference of the complex eigenvalues jUj can lead to suppression of odd versus even values. In the 
fermionic description, this is a result of BCS-type pairing in the Hamiltonian of Eq. 12.51 The anomalous term akCL-k 
and its hermitian conjugate describe pair creation and annhilation processes which might favor even occupation. We 
note that the appearance of independent Bernoulli variables for the fermion number distribution has been studied 
in the mathematical literature (see Ref. [t^) as determinantal point processes. However, it appears the cases that 
have been studied correspond to states with no BCS-type pairing and purely real eigenvalues in. These heuristic 
considerations take into account only one or two fermionic levels. Determining the global structure of P„(m) is a 
many-body problem involving a large number of fermionic levels which we address in the next section. 



B. Global structure 



Before discussing the global structure of P{Mz), we first argue that it gives a measure of the amount of disorder 
in the system. Physically this is not surprising since Sz{i) tends to disorder spins along the x direction by inducing 
tunneling between the two polarizations. This competes with the Ising term Sx{i)Sx{i + 1) which favors ordering of 
spins along the x direction. A duality transformation [t^ makes this connection more concrete. This transformation 
maps the spin Hamiltonain of Eq. 12.11 to a spin Hamiltonian of the same form but with the coupling constant 
exchanged 

N 

H = -J^ [25T,(i)T,(z + + r,(i)] . (2.25) 
1=1 

The Ta{i) operators describe dual spins that obey the same commutation relations as Sa{i) and are given by 

T,{i) = 2Sx{i-l)Sx(i) (2.26) 
T.{^ = \\\['iS.{z)]. (2.27) 
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m 



FIG. 1: Distribution function Pn(m) for g = 2, g = 0.5 and n — 50. Solid lines denote the small g binomial distribution and 
dashed lines denote the large g poissonian distribution. Notice the separate distributions for even and odd m for g > Qc- The 
center and width of the peak tend to zero for g ^ oo but approach a constant for g 0. 



in terms of the Sa{i) operators. Under this transformation the operator pn maps to 

'^«=E(^-2r,(z-l)T,(*)^ (2.28) 

where the term under the summation has eigenvalue for two dual spins ahgned on the x direction and eigenvalue 
for anti-ahgned. This implies that p„ is the number operator for dual domain walls and P„(to) is its probability 
distribution. 

In particular, the fermionic operators Ci (c|) create (annhilate) dual domain walls. For the BCS-type Hamiltonian 
of Eq. 12.51 Jg acts as the chemical potential for free fermions with dispersion — Jcos(fc) with the anomalous terms 
describing pair creation and annhilation. When g is small, the chemical potential crosses the band and pair fluctuations 
cost little energy leading to dual domain wall proliferation. When g is large, the chemical potential is below the bottom 
of the band and dual domain walls are confined due to energy cost. Note this is the exact opposite of the physical 
domain walls for the Sx{i) spins but the underlying physics of confinement versus proliferation [74] is the same. 

Analysis of Pnim) shows how the distribution function reflects the above physics Numerically, Eq. 12.171 allows for 
efficient calculation of the distribution function p„(m). The determinant representation gives XnW as a polynomial of 
degree n in e^^ and the coefficient of e™^ can be extracted via a computer algebra system to give P„ (m) . Analytical 
insight is also possible in the large-n limit through the use of Toeplitz determinant asymptotics outlined in Appendix 
IXI We find that the generating function has the following asymptotic behavior 

X„(A) ~exp[co(e*^)n + ci(e'^)logn + C2(e*^)] (2.29) 

where the coefficients Ci{z) are explicitly known functions. 

We first consider the small and large g regimes. We carry out the expansions in Appendix ^ and find for the 
generating function 

^"(^^'~^]expfe(e2^^-l) + ^(e»^-l)l - - ^^.SO) 



gc 



where we recognize the form of the generating functions as that of a binomial random variable for g gc and a sum 
of two Poissonian random variables for g 3> gc- However, notice in the g ^ gc regime that one poissonian variable 
depends on e^*^ whereas the other depends on e'^. This implies the operator pn behaves as 

'-|2P(5^)+P(^) 

where B{n,p) is a binomial random variable with n trials and probability of success p and P{X) is a Poissonian random 
variable with parameter A. Explicit expressions for P„(m) for these two regimes are given in Appendix^ We plot 
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(a) 



(b) 



FIG. 2: Calculating P{Mz) for a block of spins (red rectangle) maps to counting dual domain walls (blue circles) within. In 
the proliferated phase for g < gc (a), individual domain walls are free whereas in the confined phase (b), pairs of domain 
walls are free giving qualitatively different statistics. Notice pairs in the bulk and at the boundary contribute differently in the 
confined phase. The small arrows indicate the corresponding dual spin configuration. 



the distribution P„(m) for g = 0.5 and g = 2 in Fig. ^ Notice the distinct sphtting between even and odd values of 
m for g < gc that is absent for g > gc- 

The qualitatively different behavior has a simple physical interpretation in terms of proliferation versus confinement 
of dual domain walls. We have argued Mz simply counts the number of dual domain walls. When they proliferate 
for small g, individual domain walls are free. Each site contributes independently with probability Pqcc = 1/2 — g/A 
of being occupied and Pempty = 1 — Pocc of being empty. This leads to the binomial counting statistics that does not 
distinguish between even and odd total occupation. In contrast, the confined regime for g > gc has dual domain walls 
tightly bound into pairs that occur with low probability Ppair = 1/16(?^. In the dilute limit, the counting statistics 
of pairs is Poissonian. However, pairs in the bulk contribute a factor of two to the number of dual domain walls with 
a probability {n — l)Ppair that scales with n whereas pairs at the boundary contribute a factor of one to the total 
number with a probability 2Ppair since there are two boundaries. The differing contriubtions of pairs in the bulk and 
boundary leads to the splitting between even and odd total occupation. A schematic illustration of counting dual 
domain walls in these two regimes is illustrated in Fig. |21 

The above analysis relies on an expansion in the small parameter g or g^^ and identifying the resulting form of 
Xn(A). In order to study the evolution from the small to large g regimes, we use a large-n saddle-point integration for 
Pn{m) which works for all g. This comes at the expense of obtaining only the asymptotic gaussian behavior and is 
carried out in Appendix [CI The gaussian nature of the resulting distribution can be established two ways. We have 
already argued that p„ is the sum of n independent Bernoulli-like random variables with finite variance which means 
Pn{m) is guassian by the central limit theorem. In addition, we can examine the Taylor series 

oo . p 

\ogXn{\) =Y.-{pl)c. (2.32) 

p=0 ' 

whose coefficients (pj^)c give the cumulants or connected moments. We find all cumulants arc finite and proportional 
to n at leading order. By considering the ratio 



{pi 



\P/2 



-p/2 



(2.33) 



we conclude that the first two cumulants completely characterize the distribution for large-n which implies that P„ (m) 
is gaussian. This is not in contradiction to the small and large g regimes since the resulting distributions are also 
gaussian for sufhciently large n. 

The absence or presence of an even versus odd splitting characterizes the small and large g regimes, respectively. 
This aspect of the global structure of P„(m) reflects higher order moments which encode the underlying correlations 
of the model and appear in the large-n gaussian distributions. In particular, we find 

^-i..-n,f/2na^ [ 1 + (-1)™(1 - 5-^)^/^ 9 > 9c 

Pnim) { 1 + (-l)'"C7^-i/4 g ^ (2.34) 

[1 9 <9c 

which is of the form of either a single gaussian or the sum of two gaussians. The (—1)™ term gives the even versus odd 
splitting which has a sharp transition and critical behavior a,t g = g^ This illustrates clearly that the quantum phase 
transition manifests itself in higher order correlations of and is accessible via the distribution function P{Mz). 
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Having disucssed the physical origin of the this sphtting, we now show how it also reflects correlations of the order 
parameter or longitudinal spins Sx{i)- At first this is suprising since Pn{m) describes correlations of the transverse 
spins Sz{i)- However, recall the Sz{i) correspond to domain walls for Tx(i) which are related to Sx(i) via duality. 
Explicitly, we obtain From Eq. 12.281 

exp(i7rp„) = 4r,(0)r,(n). (2.35) 

This gives the two-point function of the dual spins Tx(0)Tx{n) as a function of pn implying that Pn{m) can be used 
to compute its average in the ground state 

(T.(O)r.(n)) = i^P„(™)(-l)™. (2.36) 

m 

The interpreation of this equation is also illustrated in Fig. 12 An even (odd) number of dual domain walls between 
two sites gives a spin configuration with spins aligned (anti-aligned) along the x direction. By using the gaussian form 
of Pn{m) in Eq. 12.341 and passing from summation to integration we find 

(r,(0)T,(n)) ~ ^ { Cn-^l^ g = (2.37) 



where we have dropped the contribution from terms entering with (—1)™ in the summation. Here C — 2^^^'^e^^^A~^^ 
with A being Glaisher's constant. This agrees with known results on the spin correlation function Sx{0)Sx{n) 
after applying a duality transformation g — > g"^. Why the agreement is exact is discussed in Appendix ICl 

We now discuss the mean fj, and variance of the guassian distributions in Eq. 12.341 They are expressed in terms 
of the Bogoliubov angle by 

We plot the mean /i and standard deviation a in Fig. |31 As 5 ^ 0, the distribution has a finite width a and a finite 
mean /i. In contrast, as g — > 00, both a and fi go to zero. This behavior can also be seen in Fig. ^ We plot Pn(jn) in 
rescaled variables in Fig. 01 The asymptotic convergence to gaussian distributions is quite good even at the critical 
point g — gc for sufficiently large n. In addition, the distinctive even versus odd splitting is absent for g < g,., present 
for g > gc, and decays algebraically for g — gc- 

Here we compare our results with that of Ref. [fi^ . Our new contribution is the demonstration and analysis of 
the even versus odd splitting in P„(m). Fig. |21 shows that this is fundamentally a boundary effect of a local block of 
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(a) 



(b) 



(c) 

FIG. 4: Distribution function Pn{m) for g = 1.25 (a), g — gc = 1 (b), and g = 0.8 (c) plotted as a function of the rescaled 
variables x = (m — {pn) c) / {p^)V^ and P{x)dx = Pn(m){p'^)l^'^ where {. . .)c denotes the cumulant. The solid lines indicate the 
asymptotic gaussian form of p„(m). Notice the convergence to guassian behavior even a.t g = g^ and the splitting between even 
and odd m for g > gc- 



spins. Thus it is necessary to separate the system size N from the block size n to see this effect. In contrast, Ref. 
[66] implicitly set A'^ = n in the beginning of their analysis. At a more technical level, we find an approach based on 
mapping the problem to Toeplitz determinants allows for separation of the length scales N and n. This allows us to 
obtain the first finite size correction to the leading gaussian behavior to P„ (m) as discussed in Appendix ^ which 
turns out to be crucial in obtaining the splitting. In addition, this approach also facilitates study of the small and 
large g regimes where the underlying physics is most evident. As another technical note, Ref. |63 also obtains the 
asymptotic gaussian behavior via a saddle-point analysis, but we show that there is an additional saddle-point which 
gives rise to the splitting and is discussed in Appendix [CI Having discussed one application of quantum noise analysis 
in a specific spin system, we now turn our attention to a discussion of general issues in experimental realizations. 
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III. EXPERIMENTAL ISSUES 



Although the use of ultracold atoms in experimental studies of quantum many-body systems has seen great progress, 
techniques based on quantum noise analysis are just becoming within reach. The fundamental difhculties arise from 
the need to collect accurate statistics on noise fluctuations of an observable, ideally to high orders. Here we discuss 
some proposals for experiments and the relevant issues for quantum noise analysis. 

Optical lattices offer one of the most promising avenues for engineering many-body systems in ultracold atoms 
including spin systems [85l | but they may pose problems for noise analysis. This is especially true for low-dimensional 
systems. They have already been realized in optical lattices |2l|, |22, |5g| but occur in a large array of up to thousands 
of one- or two-dimensional systems within each sample. A single measurement such as an absorption image over 
the entire sample corresponds essentially to ave ragi ng over the large array which supresses fluctuations. Low-order 
statistics such as second-order coherence |H3)IE3jlH3 still accessible. However higher order noise fluctuations may 
be difficult to detect as a result of the intrinsic averaging that occurs within each sample. 

Accessing higher order noise may be possible through realizations of single or small numbers of low-dimensional 
systems where fluctuations are not negligble compared to the average. Such an approach has already been used to 
study the Berezinskii-Kosterlitz-Thouless transition |63j . In this experiment, the quantum noise in the interference 
amplitude of just a pair of two-dimensional quasi-condensates was used to detect signatures of the transition. For 
one-dimensional systems, surface microtraps allow for tight confinement in a wa veg uide geometry |75l l76l and 
optical box traps have recently demonstrated single one-dimensional condensates |78l |. 

In addition to preparation, experiments will also have to advance in measurement techniques. The ideal measure- 
ment of quantum noise for magnetization would be state-selective in situ imaging with single atom precision to obtain 
accurate high order statistics and high resolution to resolve spatial correlations. This has not yet been achieved with 
optical absorption imaging but may be possible through proposed scanning tunneling_ microscropy of ultracold atoms 
|79j. However, single atom counting has been demonstrated in an atom laser setup ^Ij. However, this is more suited 
to the study of temporal correlations rather than the spatial correlations describing noise in magnetization. Detection 
of single atoms has also been used in Hanbury-Brown-Twiss experiments but with imaging after expansion which 
resolves momentum correlations and not in situ which probes real space correlations. An alternative measurement 
involves imprinting the quantum noise for magnetization onto a weakly coupled laser probe 80]. In this case, the 
quantum noise of the laser itself should be less than that of the magnetization and may make features at the single 
atom level difficult to detect. 

The above discussion of general issues for quantum noise experiments also applies to the transverse Ising model 
we considered. However, resolving the even versus odd splitting in the distribution function requires additional 
considerations. As shown in Fig. [21 distinguishing between even and odd magnetization requires a sharp boundary 
between the block of spins undergoing measurement and the rest of the spin chain. This requirement does not pose 
difficulties for measurement via high resolution in situ imaging. However, for techniques that couple to the entire 
sample, RF cutting (see Ref. 81J) may be needed to isolate the region of interest beforehand. 

Sensitivity at the single atom level will also be crucial for detecting the even versus odd splitting. We consider 
the feasibility of such a measurement in the absence of in situ imaging by analyzing the proposal of Ref. [^O] for 
coupling the magnetization to off-resonant coherent laser light. For our purposes, the main result relates the output 
quadrature of light Xout to the input quadrature Xin and the magnetization as 

Xout ^ X,r. - (3.1) 



where k is the atom-photon interaction and na is the number of atoms. Here the input quadrature component XiU 
has been rescaled by the number of photons Up and describes coherent light with noise AX in = l/\/2 

The quantum noise of the magnetization is thus imprinted on top of the inherent photon shot noise of light. 
Discerning gross features of the magnetization distribution only requires fluctuations in Xout due to quantum noise in 
the collective magnetization be greater than fluctuations due to photon shot noise. The width of the magnetization 
distribution gives an estimate AM^ ^ y^n^cr where a is shown in Fig. |31 and yields the condition 

This enables resolving the smooth distribution shown in Fig. ^ iu the regime with no splitting whereas it only gives 
the average of the two separate distributions of the same flgure in the regime with splitting. Detecting this splitting 
requires Xout fluctuations due to single spins rise above photon shot noise and yields a much more stringent condition 

^>M- (3.3) 
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Eq. 13. II describing the coupling of magnetization to light is valid when decoherence due to spontaneous emission is 
small and for low photon absorption. This implies the atomic depumping r] — UparT^ / AA^ and photon absorption 
e = na<Jr^^ / AA^ are small. Here Up is the number of photons, ct^ is the resonant scattering cross section which is 
distint from a the width of the magnetization distribution, T is the spontaneous decay rate, A is the illuminated 
cross sectional area, and A is the detuning. The atom-photon interaction is given by n = ^yn^n^ar^ / AA. For given 
•q,e <^ 1, the condition for resolving collective magnetization noise in Eq. 13.21 becomes 



- < aV2^ (3.4) 
which can be satisfied for sufficiently large detuning. However, Eq. 13.31 for resolving single spin noise becomes 

which is likely an unphysical regime since since Ur ~ where A is the laser wavelength and is a small parameter. 
This states that the photon shot noise for coherent light overwhelms the single spin noise. Using squeezed states of 
light with smaller quadrature fluctuations may allow detection of single spin noise but a more detailed analysis would 
be necessary. 



IV. CONCLUSION 



Ultracold atoms offer the prospect of experimentally realizing quantum many-body systems in a strongly interacting 
regime with advancements such as optical lattices, Feshbach resonances, and confinement to low dimensionality. Noise 
analysis, or the study of fluctuations of an operator potentially to high orders instead of just the average, offers an 
avenue to access correlation functions and to characterize many-body states. We discussed a proposal for using noise 
analysis in quantum spin systems realized with ultracold atoms. By studying the entire distribution function for the 
magnetization of a mesoscopic spin system, the behavior of high order correlation functions can be extracted. We 
analyzed the Ising spin chain in a transverse field in detail, focusing on the transverse magnetization. The resulting 
distribution function possesses a distinct splitting between even and odd values of the magnetization. This splitting 
behaves differently in ordered, critical, and disorderd phases and has a simple interpretation in terms of confinement 
versus proliferation of domain walls. Finally, we discussed pertinent issues for implementing quantum noise analysis 
through experiments in the near future. 
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APPENDIX A: TOEPLITZ DETERMINANTS 



In this appendix, we outline the derivation of the asymptotic form of the generating function in Eq. 12.291 using 
the theory of Toeplitz determinants. Eq. 12.171 gives an exact expression for the Xn{m) as a determinant of a matrix 
with constant diagonals. The matrix elements are given by the fourier coefficients fi_j(z) of the function f{z,x) in 
Eq. 12.181 The theory of Toeplitz determinants (see Ref. j82|) gives the large-n asymptotics in the form of 12.291 with 
applications to a number of physical problems (see Ref. [83 and references therein). 

The simplest case to consider is when f(z,x) is continuous on \x\ = 1 with integer winding number where a variant 
of Szego's lemma holds For power-law singularities in f (z. x) (branch cuts in log f{z,x)), the Fisher-Hartwig 

conjecture is needed and has been proven in the cases we need |84| . These results are rigorous but we can give a 
heuristic motivation for their structure. Formally, one can write 

det[/,_,(z)] = exp(Tr[log/,_,(z)]) (A.l) 

where i,j = 1 . . .n and the leading contribution comes from the replacement 

lik 

^-log/(.,e*'=). (A.2) 
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This replacement corresponds to noticing fi^j{z) is translationally invariant in the sense that it depends only oni — j 
and evaluating the trace in momentum space. The subleading terms correspond to finite size corrections in n. 

The functions Ci{z) depend on the analytic structure of log/(2;,a;) and differ for g < Qci 9 = So and g > gc- It will 
be useful to define the following functional acting on functions of x 

holfix)] = ^/^log/(^) (A-3) 

Kim] = 

where the contour integrals are along the unit circle and notice hy[f{x)] depends on an additional parameter y. If 
there are singularities on the unit circle then the integrals are in the principal value sense. 

First we consider g > gc- For all z, f{z, x) is continuous on \x\ = 1 with zero winding number. We obtain 

co{z) = ho[f{z,x)] (A.6) 
ci{z) = (A.7) 
C2{z) = E[f{z,x)]. (A.8) 

Next we consider g < gc- For Re[z] > (Re [2] < 0), log/(^,a;) is non-singular with winding number zero (one). 
We obtain 

r (z^ = J'^ot/^^'^)] Re[z]>0 
^ \ho[x-'f{z,x)]+\og{-0 Re[z]<0 
ci{z) = (A.IO) 

^ \E[x-'f{z,x)]+2h^[x-\f{z,x)]-\og§^{z,0 Re[z]<0 ' 
where |^| < 1 corresponds to the singularity of log f(z,x) closest to the unit circle given by 



^_ 2ff + V2-z^-z-2-V 

and the square root is positive for positive real arugments. 

The final case we consider is g = gc- Here, f{z,x) is discontinuous at a; = — 1 for all z. This singularity can be 
written in power-law form 



where the exponent is given by 

iy(z\ 



a{z) = 7 - - arctan(2). (A. 14) 



Notice log/(^;,a;)/gf(^,a;) is continuous. We obtain 

co{z) = ho[f{z,x)/g{z,x)] (A.15) 
ci{z) = -aizf (A.16) 
C2{z) = E[f{z,x)/g{z,x)] + 2a{z)h.i[f{z,x)/g{z,x)]+\og[Gil+aiz))Gil-a{z))] (A.17) 

where G{x) is the Barnes-G function satisfying G{x + 1) = r(a;)G(a; -|- 1) with r(a;) the gamma function. 
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APPENDIX B: SMALL AND LARGE g LIMITS 



We now turn to the small and large g expansions of the generating function XnW in Eq. 12.301 For g ^ gc, we use 
the determinant representation in Eq. 12. 171 and first expand the matrix elements fn{z) in powers of g. We find 



fn(z) 



9^ + 0{g') 
-g^ + 0{g^) 
0{g') 



n = 1 
n = 2 
n 7^ 0,1,2 



(B.l) 



The determinant det[/i_j(z)] is lower triangular and given by the product of the diagonal matrix elements 



XniX) 



2 + .g ^ 



(B.2) 



as in Eq. 12.301 This corresponds to the generating function for a binomial random variable and the identification 
of Pn ~ B{n, (2 — 5)/4) in Eq. 12.311 A binomial variable B{n,p) with n trials and probability of success p has a 
distribution 



Prob[i?(rt,p) = m] = 



r(n + i) 



r(m + l)r(n-m + l 



■(1 



(B.3) 



where Tin) is the gamma function. The distribution P„(to) is given by Eq. IB.3l with p = (2 — .g)/4 which we use to 
plot Fig. Ill ' 

For g ^ gc, we analyze the determinant of Eq. 12.301 using the large n asymptotic of Eq. 12.291 as discussed in 
Appendix ^ Instead of the integral representation employed in Eq. IA.6I it will be easier to use the equivalent 
representation 



ho{z)n + ^ khk{z)h_k{z) 



k=l 



X„(A) ~ exp 

in terms of the fourier coefficients of log/(z, x) given by 

oo 

^og f{z,x) = ^ h,i{z)x". 

n— — OO 

We expand the coefficients hn{z) in powers of g and find 

r^(z2-i) + o(g-3) ^ = 

hn{z)^ <Ti-g{z-l) + 0{g-^) n = ±l 

l0(<?-3) n^O,±l 



(B.4) 



(B.5) 



(B.6) 



Using Eq. IB. 51 this gives the generating function to 0{g ^) as 

n — 1 



X„(A) = exp 



I652 



tig- 



(B.7) 



which is that of the sum of two Poissonian random variables as in Eq. 12.301 However, notice one term enters with 
g2iA g^^^ ^j^g other enters with e'^ meaning we identify p„ ~ 2P((n — l)/lQg^) + P{l/8g^) in Eq. 12.311 A Poissonian 
variable P(A) with parameter A has a distribution 



Prob[P(A) = m] = 



;-^A" 



r(m + 1) 

and a sum of two Poissonian variables of the form 

X = 2P(A2) + P(Ai) 



(B., 



(B.9) 
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has a different distribution for even and odd values 

m 

Prob[X = 2m] ^ Prob[P(A2) = n] Prob[P(Ai) = 2(m - n)] (B.IO) 

n=0 

m 

Prob[X = 2to + 1] = ^Prob[P(A2) = n]Prob[P(Ai) = 2(m-n) + 1]. (B.ll) 

n=0 

These sums can be evahiated in terms of the hypergeometric function iPi(a; b; z) 

(2A2 + Ai)^] / ("^5 5; ^) 



Prob[X = m] — exp 



4A2 



A2 <( '"^ — i . ^ . _ ^ (B.12) 



The distribution P„{'m) is given by Eq. iRTl with Ai = l/8.g^ A2 = (n - 1)/I6.g2 which we use to plot Fig. □ 

APPENDIX C: SADDLE POINT ANALYSIS 

In this appendix, we obtain the distribution P„(m) by evaluating the integral of Eq. 12.81 in the saddle-point 
approximation. We rewrite the integral as a contour integral over the unit circle |z| = 1 with the asymptotic form of 
X„(A) in Eq. TTM to obtain 

1 f dz 

^" ^ f ^;^6^P['^o(z)« + ci(z)logn + C2(z)]. (C.l) 

The leading contribution to integral for large n is controlled by the term z~™ exp[co(z)n] which has saddle points at 

m = hzcq^z) (C.2) 

where primes denote differentiation with respect to z. For the function co{z) obtained in Appendix 1X1 we find there 
are two saddle points. One is located on the positive real z axis and the other is on the negative real z axis which 
we denote as z± respectively. In addition, the |z| = 1 contours can be deformed to steepest descent contours and we 
obtain in the saddle-point approximation 

, ^ 2;±^"+^^exp[co(2;±)n + ci(2±)log7i-hc2(z±)] 
± ^2T:[{m + l)/zl+nc'^{z±)] 

We have already argued that the in the large-n limit, Pnim) should be gaussian in the discussion of Eq. \2.',VA\ In 
general, the gaussian form of P„(m) is good near the peak of the distribution. Thus, we will perform an expansion 
of Eq. IC.3I about the peak located at n^± and the corresponding saddle-points at . The two quantities are related 
by the zeroth order saddle-point condition 

^l = zX{z°) (C.4) 

and we have dropped the ± subscript for clarity. We then substitute 

m = n/.t + nSm (C.5) 
z = z" + dz (C.6) 

into the full saddle-point equation of Eq. IC.2I and solve order by order to obtain 

Sz — ai6m + a2Sm^ (C.7) 



where the coefficients Ui are given by 
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Evaluating Eq. IC.3l up to second order in Sm for the leading term z ex'p[cQ{z)n] and zeroth order for the other 
terms yields 

±1 p-(™-"A'±)^/2"cr±-™log2± 

P„(m) ~ — ^= exp [co(4)n + ci{zl) logn + C2{z°±)] (C.IO) 



where CTj. is given by 



a'^z'[c'o(z^) + zX{z")] (C.ll) 



and the ± subscripts have been suppressed. Since nfi should correspond to the peak of the distribution, we see that 
is determined by the condition that the term linear in m of the exponent in Eq. IC.IOI should have a real part of 
zero to ensure that it does not shift the peak. This implies that z^. = ±1. 

At this point, we have to check whether the z^. saddle-points contribute to the leading term z~™ exp[co{z)n] 
at the same order. The functions Co{z) are given in Appendix IXI and from Eq. I A. 61 and Eq. IA.15I we find that 
Co(±l) = which means that both saddle-points should be included for g > gc- However, from Eq. IA.9I we see that 
Co(+l) > co(— 1) due to the additional contribution of log(— ^) of Eq. lA. 121 with |^| < 1. Thus the contribution from 
z° should be dropped for g < gc- 

Evaluating the remaining expressions in Eq. IC.IOI using the functions Ci{z) obtained in Appendix IXI gives the 
gaussian form of Eq. 12.341 In addition we find that n± = fJ. given by Eq. IC.4I corresponds to the mean and CTj. = 
given by Eq. IC . 1 II corresponds to the variance and are given by the expressions in Eq. 12.381 An enlightening way of 
writing Eq. IC.IOI is given by 

P„{m) 7==^ [Xn(0) + (-l)™XnW] . (C.12) 



where with A = 0,7r correspond to contributions from the zeroth order saddle points z — e'^ = +1,-1, 

respectively. The contribution from the z ~ +1 saddle point is dictated by normalization of the distribution function 

Xn(0)=^F„(TO) = l. (C.13) 

m 

However, Eq. 12.361 shows the z = — 1 contribution 

X„(7r) = ^P„(m)(-1)" = 4(r,(0)r,(n)). (C.14) 

m 

is determined by the Tx(0)Tx{n) correlation function. Since the asymptotic expression we obtained for x„(A) is exact 
for any A and the even versus odd splitting is given precisely by Xni''^), we see that the splitting gives the exact 
asymptotic behavior of {Trc{0)Trc{n)) . 
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